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ABSTRACT 

The  relaxation  time  scale  in  glassy  materials  is  derived  within  a  model  of  anomalous 
defect  diffusion.  The  effect  of  the  defects  on  ion-doped  polymeric  glasses  is  to  produce  a 
stretched  exponential  waiting  time  distribution  for  ion  jumps.  The  characteristic  time  scale  for 
ion  jumps  is  connected  the  temperature  and  pressure  dependent  concentration  of  mobile  defects. 
The  resultant  expression  for  ionic  conductivity  is  compared  with  experimental  results  for  the 
polymer  electrolytes  poly(ethylene  glycol)  (PEG)  and  poly(propylene  glycol)  (PPG)  containing 
LiCFaSO^. 


1.  INTRODUCTION 

In  1 889,  Arrhenius  introduced  the  concept  of  an  activation  energy  together  with  a  law  for 
the  related  characteristic  relaxation  time  Qxp{AlkT)  where  k  is  Boltzmann's  constant 

and  T  is  the  absolute  temperature  in  degrees  Kelvin.  The  Arrhenius  law  was  later  derived  by 
Kramers  in  terms  of  the  trajectory  of  a  particle  successfully  crossing  an  energy  barrier  of  height, 

A ,  with  an  attempt  frequency  of  .  However  the  Arrhenius  law  is  typically  not  valid  for  most 
glass-forming  materials.  A  corresponding  law  for  these  materials  was  proposed  by  Vogel^  in 
1921,  for  viscosity  77  in  the  form  //(T)  =  77  ^00  ^  Vogel  used  this  empirical  fit 

OC 


for  viscosity  experiments  on  mercury,  water,  and  oils.  In  terms  of  a  time  scale  Vogel’s  law  can 
be  written  equivalently  and  more  transparently  as, 

B.. 


r,,  =  /!„  exp 


T-T. 


(1) 


oy 


This  later  form  was  proposed  by  Fulcher^  in  1925,  and  Tammann  and  Hesse^  in  1926.  Today 
this  is  called  the  Vogel  law  or  the  VFT  law  or  some  other  combination  of  the  above.  The 
interpretation  of  the  parameters  Ay,  Bv  and  is  not  so  straightforward.  First,  one  notices  that 

Tq  is  a  special  temperature  where  the  time  scale  diverges.  The  temperature  Tq  is  typically  well 
below  'T. ,  the  glass  transition  temperature  of  the  material,  /.e.,  the  relaxation  dynamics  are 


focused  on  and  not  on  .  There  have  been  several  attempts  to  derive  the  Vogel  law  or 
alternative  laws.  Previously,  we  had  derived  the  equation  for  the  temperature  dependence  of  the 


relaxation  time  at  zero  pressure 

~  ^UD 


(T-Tj'- 


(2) 


where  A^^,Bi^,^and  are  constants.  The  derivation  of  eq.  (2)  starts  by  initially  arriving  at  a 

stretched  exponential  relaxation  time  distribution  within  a  defect  diffusion  model.  The 
relaxation  time  scale  is  the  mean  time  of  the  stretched  exponential  law.  This  time  scale  depends 
on  the  (temperature  and  pressure  dependent)  mobile  defect  concentration.  In  this  model,  the 
defects  carry  free  volume  and  can  unfreeze  the  parts  of  the  glass  that  they  visit.  With  enough 
defects  the  glass  is  not  frozen.  As  the  temperature  is  lowered,  defects  cluster,  and  the  number  of 
mobile  defects  decreases.  As  this  occurs,  the  material  becomes  more  viscous  (rigidity  begins  to 
set  in).  At  ,  the  defect  concentration  decreases  to  the  point  where  rigidity  percolates,  and  the 

glassy  state  is  formed.  Relaxation,  however,  is  still  occurring.  A  phase  transition  in  the  number 
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of  mobile  defects,  as  the  temperature  is  lowered  towards  ,  creates  the  behavior  characterized 
by  eq.  (2). 

Equation  (2)  has  consistently  been  as  good  as  or  better  than  the  Vogel  law  (eq.  1)  for 
fitting  ionic  conductivity,  dielectric  relaxation  and  viscosity  data  for  glass-forming  materials  '  . 
The  Vogel  law  often  fails  to  fit  data  adequately  near  and  an  Arrhenius  law  is  sometimes  used 

near  ,  in  conjunction  with  a  Vogel  law  in  the  region  above  .  Equation  (2)  provides  a 

consistently  better  fit  throughout  both  regions. 

In  this  paper  we  have  extended  eq.  (2)  to  include  the  effect  of  pressure.  It  will  be 
demonstrated  that  the  resultant  generalized  Vogel-type  law,  eq.  (6),  successfully  describes 
measurements  of  ionic  conductivity  in  salt-containing,  glass-forming  liquids,  such  as  those 
frequently  used  as  the  electrolyte  in  a  battery, 

II.  THEORY 

Consider  a  glass-forming  material  possessing  a  concentration  c  of  defects,  where  c,„  of 
these  are  mobile.  Ion  hopping  motion  at  an  ion  occupied  site  occurs,  at  time  t,  due  to  the  flux  of 
defects,  F(l),  into  that  site.  To  calculate  F(t)  we  employ  a  waiting  time  probability  density,  y/it), 
for  holding  a  defect  for  a  time  /,  in  between  its  Jumps.  This  incorporates  a  randomness  into  the 
defect  motion.  The  relaxation  function  ^(/)  is  given  by  =  exp(-F(/)) .  The  defect  flux  into 
a  site,  within  a  time  t,  is  c„J^(t)  where  N(t)  is  the  number  of  distinct  sites  a  random  walking 
defect  visits  within  a  time  t.  The  relaxation  law  is 

$>(/)=  exp(-c„,V(/)).  (3) 

For  a  typical  Brownian  motion  type  random  walk,  N(t)  is  proportional  to  /  and  exponential 
relaxation  occurs.  Slower  than  exponential  relaxation,  i.e.  stretched  exponential  relaxation  is 
another  possible  relaxation  law,  i.e. 
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(4) 


(t){t)  =  exp(-  Ac^  t")  =  exp(-  [/  /  rf ) 
where  (3  <  \  and  /I  is  a  constant  and  on  the  RHS  r  =  .  This  is  the  typical 

behavior  exhibited  by  glassy  materials.  There  are  several  ways  in  which  one  can  arrive  at  eq. 
(4).  If  one  relates  a  time  scale  to  overcoming  a  free  energy  barrier  as  follows: 
/  = /q  exp(- ( A  -  / ^7’)  then  the  distributions  of  energy  barriers,  A,  entropies,  S,  and 

prefactors,  to,  can  each  generate  a  distribution  y/{t)  of  waiting  times  between  jumps.  Consider 
first  that  only  S  is  a  random  variable  and  the  distribution  of  environments  produces  is 
/(5)  =  exp(-6'/5'o) .  Then  =  f{S)dS  leads  at  long  times  to  l/r'”''’ with 

{3  =  k/S^.  When  S<Sq,  /?  <  1 .  For  this  case,  for  a  single  defect,  the  mean  waiting  time 

CO 

between  Jumps,  is  infinite.  If,  instead,  one  places  all  the  randomness  in  the  energy 

(.) 

term  a  similar  stretched  exponential  decay  arises,  but  with  the  exponent  being  temperature 
dependent.  In  general,  both  mechanisms  can  be  expected. 

Since  only  the  defects  are  mobile  in  the  model,  as  the  temperature  is  lowered,  the  defects 
cluster  (or  correlate  their  motion)  to  lower  the  system  entropy.  We  now  make  the  assumption 
that  single  defects,  of  concentration  c, ,  are  more  mobile  than  a  cluster  of  defects.  We  therefore 
replace  c,„  in  eq.  (4)  by  c/.  To  have  a  single  (isolated)  defect  at  a  site,  one  must  first  have  a 
defect  there  with  probability  c  and  also  have  all  of  the  z  neighbor  sites  within  its  correlation 
volume  unoccupied,  i.e. 

c,  =  c(l  -  c)'  (5) 

with  z  =  {i^/dy ,  where  ^  is  the  defect-defect  correlation  length,  and  d  is  the  nearest-neighbor 

lattice  spacing.  In  a  mean  field  lattice  gas  model,  the  correlation  length  ^  between  the  defects 
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grows  near  and  above  a  critical  temperature  T^as  ^{T)k  L 


where  L  is  a  constant  and 


Tc  is  the  temperature  at  which  single  defects  disappear  and  below  which  only  defect  clusters 
exist.  With  increasing  pressure,  the  nearest-neighbor  spacing  d  is  assumed  to  decrease 

isotropically  as  ^  d^\\-6{T,P))  where  \-5{T,P)  =  V{T,P)IV{T,Qi)  is  the  fractional 

volume  change  of  the  material  as  pressure  increases  and  do  is  the  mean  lattice  spacing  separation 
at  zero  pressure.  The  time  scale  in  the  stretched  exponential  can  now  be  expressed  as, 


cr'''^ro  =c- 


Tf,  exp 


{T-Ty'-{\~5{T^P)) 


(6) 


where  B  =  -(L  do)''  \r\(I-c)  f>.  This  represents  a  new  relaxation  law  that  is  Vogel-like,  but  with  a 
3/2  temperature  exponent  and  the  inclusion  of  pressure  effects.  Note  that  Tc  is  a  function  o^P. 

In  an  ion-containing  polymeric  glass-forming  material,  as  described  by  our  model,  ion 
transport  is  controlled  by  the  defects  and  it  is  assumed  that  the  relaxation  times  for  ionic 
conductivity  are  the  same  as  those  for  dielectric  relaxation.  Consequently,  a  model  of 
conductivity  can  be  developed  as  follows.  Although  an  individual  defect  has  an  infinite  waiting 
time  between  jumps,  the  ion,  hit  by  a  flux  of  defects,  has  the  stretched  exponential  waiting  time 
distribution.  All  the  temporal  moments  of  the  stretched  exponential  are  finite,  including  the  first 

moment,  which  we  label  as  xi.  Since,  the  diffusion  constant  is  of  the  formD  =  - — ,  employing 
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Basically,  as  the  pressure  is  increased,  defects  are  pushed  closer  together  and  become  more 
clustered  leaving  fewer  single  (mobile)  defects.  This  will  decrease  the  defect  flux,  increasing  the 
time  scale  x,  and  thus  decreasing  the  conductivity  a. 

An  equation  for  the  viscosity,  tj,  follows  by  applying  the  Stokes-Einstein 
kT 

equation  rj  - - ,  where  ro  is  the  effective  ion  radius. 


III.  COMPARISON  WITH  EXPERIMENT 

For  the  purpose  of  comparing  the  theory  of  the  ionic  conductivity  with  experiment,  eq. 
(7)  is  written  as  follows: 


(y{T,P,)  = 


T{\-S) 


exp 


BT 


{T-TJ'\\-5) 


(8) 


where  and  B  are  constants.  The  S  in  the  denominator  of  the  pre-exponential  accounts  for  the 

increase  in  ion  concentration  as  pressure  increases.  Other  terms  in  the  pre-exponential  may  be 
pressure  dependent  but  are  taken  to  be  constant  in  the  present  paper.  It  is  assumed  that  the 
volume  of  the  material  changes  with  pressure  according  to; 

1-^  =  1-  x{T)P  +  f{T)P^  +  g{T)P^  ■  (9) 

This  form  is  chosen  as  it  is  supported  by  the  PVT  data  of  Zoller  and  Walsh.  Finally,  the 
critical  temperature  is  assumed  to  be  pressure  dependent . 


To  summarize  the  results  for  ionic  conductivity  vs.  pressure,  the  theoretical  values  for  In  a 
vs.  pressure  at  three  temperatures  were  calculated  using  eq.  (8).  The  experimental  values  were 
calculated  using  eq.  (12).  The  theoretical  and  experimental  values  are  shown  in  Fig.  1.  The 
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Figure  shows  good  agreement  between  theory  and  experiment.  The  investigation  of  the  first  and 
second  pressure  derivatives  of  the  conductivity  will  be  reported  elsewhere. 

IV.  CONCLUSIONS 

A  defect  diffusion  model  has  been  developed  that  is  capable  of  describing  the 
temperature  and  pressure  variation  of  several  dynamical  processes  in  glass-forming  materials. 
The  basic  principle  of  the  model  is  that  electrical  relaxation  or  ionic  conduction  occurs  when  a 
defect  encounters  a  dipole  or  ion.  Each  of  the  parameters  in  the  theory  has  a  clear,  physical 
interpretation.  There  is  an  underlying  temperature  that  represents  the  temperature  below  which 
no  defects  are  mobile.  There  is  a  dimensionless  constant  in  the  exponent  that  depends  upon  the 
correlation  length,  separation  and  concentration  of  the  defects.  For  dielectric  relaxation,  the  pre¬ 
exponential  is  a  product  of  a  characteristic  relaxation  time  (taken  to  be  about  that  for  a  lattice 
vibration)  and  the  defect  concentration.  The  pre-exponentials  for  ionic  conduction  and  viscosity 
contain  the  expected  additional  terms.  Comparison  of  theory  and  experiment  reveal  an  excellent 
representation  of  the  data  over  a  wide  range  of  temperatures  and  pressures. 
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FIGURE  CAPTIONS 

Figure  1.  Electrical  conductivity  vs.  pressure  for  PEGTiCFsSOs  and  PPGiLiCFsSOa  vs. 
temperature.  The  points  are  experimental  and  the  lines  are  theoretical. 


